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A general technique of exact calculation of any correlation functions for the special class of one- 
dimensional spin models containing small clusters of quantum spins assembled to a chain by alter- 
nating with the single Ising spins is proposed. The technique is a natural generalization of that 
in the models solved by a classical transfer matrices. The general expressions for corresponding 
matrix operators which are the key components of the technique are obtained. Despite the spatial 
inhomogeneity, such class of spin chains is proven to have only one correlation length for all types 
of correlation functions. The technique is illustrated on two examples: symmetric diamond chain 
and asymmetric sawtooth chain. 
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I. INTRODUCTION 

Recently much attention has been paid to the exactly 
solvable models of one-dimensional magnetism, consist- 
ing of the small quantum spin clusters connected to 
each other via Ising spins^— For the sake of simplicity 
we will call them Heisenberg-Ising chains (HIC). The 
interest toward such kind of spin systems stems out from 
the several features they possess. First of all, they allow 
one to obtain an exact statistical-mechanical solution 
of the problem of strongly interacting spins in terms of 
classical transfer matrbs^Zr— . In addition to that, such 
systems combine the quantum and classical properties 
and also, to some extent, are connected to the problems 
and models of molecular magnetism 3 -!, as the classical 
transfer matrix is constructed by the diagonalization 
of the small quantum spin cluster. And finally, very 
recently the family of trimetallic coordination polymer 
compounds has been sythesized 3 - 2 ^— . These compounds 
in the magnetic sense are indeed a one-dimensional 
magnets with both classical (Ising) and quantum 
(Heisen berg ) bonds. The single-chain magnet reported 
in Ref. l32| is an example of a Heisenberg-Ising chain 
of triangles with three-spin linear quantum clusters; 
another interesting molecular magnet (spin ring) based 
on the highly anisotropic properties of the Dy 3+ ions has 
been reported in Ref. l33l . Both systems have been ana- 
lyzed by the transfer-matrix solution of Heisenberg-Ising 
spin chains in Ref. H3- Another class of examples is given 
by chain magnets with two alternating magnetic ions, 
thus giving rise to alternating quantum and classical 
spins (see, e.g., Refs. [35M37b . Let us describe the general 
scheme of the HIC: 



• There is a small cluster of a few quantum spins 
with a certain topology of bonds which could 
contain various types of spin-spin interaction as 
well as interactions with external fields. The 



Hamiltonian of the single small cluster is, thus, a 
finite dimensional matrix. 



• These clusters are assembled to the chain or any 
other one-dimensional structure by alternating 
with the single Ising spins. 



• The interaction between the single Ising spins and 
any quantum spin from small cluster must contain 
only z-components (Ising interaction). 

For the mixed quantum-classical spin system with the 
above mentioned properties the exact solution is obtained 
in the following way. The idea is to apply a classical 
transfer-matrix formalism in which the entries of classical 
transfer matrix is obtained from the diagonalization of a 
small quantum spin cluster. Although the number of pa- 
per devoted to the exact solutions of the HIC and to their 
magneto-thermal properties is rather impressive, only 
few models have been solved within the classical trans- 
fer matrix formalism^ ^ 11 ' 13-17 ! 19-22 ] 26 , while the other 
method of so-called decoration iteration transformation 
has been also used!"— & 12 i 24 i 25 . Particularly, within the 
direct classical transfer matrix technique the following 
models have been solved exactly: the double tetrahedra 
chain with quantum triangles of spins&ii, the sawtooth 
chain with asymmetric pair of quantum spin o 10 ' 13 , di- 
amond chai n 14 i 16 ' 17 i 20 ~ 22 i 26 and orthogonal-dimer chain 
with anisotropic and asymmetric quantum triangles^. 
Despite the complete and exact description of the ground 
states and ground states phase diagrams, all thermody- 
namic functions, finite-temperature entropy and magne- 
tization, magnctocaloric effect, etc, the expressions for 
the correlation function are still missing. In the present 
paper we are going to fill this gap by developing the exact 
technique for the calculation of any finite-temperature 
correlation function for any HIC. 
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The paper is organized as follows. In the second Sec- 
tion we remind the reader the main features of the clas- 
sical transfer-matrix formalism for the HIC, in the next 
section we present the formalism of calculating the corre- 
lation function of HIC in terms of traces of certain matrix 
product and proof that for all kinds of correlation func- 
tions the correlation length is the same. In the Section 
IV we apply the general formalism for several models. 
The paper ends with a Conclusion. 
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(2.1) 



here a n is the Ising spin situated between n-th and (n+1)- 
th small cluster and H is the magnetic field acting on 
Ising spins. Due to the properties listed above the block 
Hamiltonians commute to each other, 



II. CLASSICAL TRANSFER-MATRIX 
FORMALISM FOR HIC 

The Hamiltonian of a HIC has the form of a sum of 
a local block-Hamiltonians H n which correspond to the 
n-th small quantum cluster and its interaction with left 



[Hi,Hj] = 0. 



(2.2) 



This fact allows one to expand the exponential in the 
expression for the partition function of the model and 
obtain an expression which by its structure corresponds 
to the classical transfer matrix scheme: 



J 



N 



Z = Tr s ^^exp(-^)=^n e "^ f(tT " +CT " +l) E' n: » e "^" 

(<r) (t) (a) n=l (r) 



N 



N 



(2.3) 



(a) n=l ( T ) t=l (<t) n=l 



here the cyclic boundary conditions are supposed, 
Xi(a n , cr n +i \T n .b) are the eigenvalues of the block Hamil- 
tonian H n (which is supposed to be an L by L matrix ) 
depending on the values of classical spins r and classi- 
cal spins a n which alternate with the blocks, Tr„ means 
trace over the states of the n-th quantum block, and Tr 
stands for the trace for the 2 by 2 transfer matrix with 
entries T a an+1 . After summation over the values of r 
for each block one gets the expression depending only on 
a n and cr ra +i which is the classical 2 by 2 transfer matrix 



(2.4) 



where Z an>!Tn+1 = £ (r) ^ti e~^">^ T ^ is the 
partition function for one block, depending on the values 
of two neighbor Ising spins, signs "+" and "— " in in- 
dexes correspond to a = 1/2 and a = —1/2 respectively. 
The formalism of classical transfer-matrix allows one to 
obtain exact analytic expressions for almost any expecta- 
tion values of the function of local variables (spins) and, 
of course, for any thermodynamic functions, as the ex- 
pression for the free energy per block can be obtained in 
a straightforward way by finding the maximal eigenvalue 
of the transfer matrix (12.41): 



7 io 4' 



PH y 



,-BH 



(2.5) 



From this expression any thermodynamic function, like 
magnetization, entropy, specific head e.t.c. can be ob- 
tained exactly. However, to obtain correlation functions 
of the spins demands more efforts and further develop- 
ment of the classical transfer- matrix technique. 



III. CORRELATION FUNCTIONS AND 
CORRELATION LENGTH 



The block Hamiltonian can contain various terms with 
quantum and classical spin variables and which are de- 
termined by the topology of the elementary block and 
the character of interactions. Let us denote the quantum 
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spins by S; then one can suppose that block Hamiltonian 
depends on S, a and may be other classical spins 77,, 



T-L n — T-L ( T n ,b) 



(3-1) 



Let us concentrate on the scheme of calculation of the ex- 
pectation values of combinations of spins operators. First 
of all, let us start form the simplest case, the expectation 
values of er-s, which is quite similar to that of the Ising 
chain^i. Thus, for the thermodynamic average of just 
one <7 one obtains 

1 N 

(a) n=l 

= -7j y ^^ T ai ^ 2 ...i7,jT C y j , crj+1 ...T ( y NtC j l . (3.2) 

This expression is the trace of the product of N transfer- 
matrices with the matrix 



1/2 
-1/2 



(3.3) 



inserted before the j-th transfer matrix in the trace for- 
mula. Thus, we have 

fo) = |Tr (a z T N ) = iTr (a z AA w A- 1 ) , (3.4) 

where we used the matrix identity T = AAA -1 , with the 
diagonal matrix A and the orthogonal matrix A. Thus, 



A = 



Ai 
A 2 I ' 



A=| COS f ~ sin / ], (3.5) 
Slllffl cos© ' 



where Ai^ are the eigenvalues of the transfer-matrix, and 
</> is a parameter. Having all these notations we now 
can express the statistical mechanical average from Eq. 
(|3.2[) in terms of the entries of the orthogonal matrix 



diagonalizing the transfer-matrix. Using the properties 
of the trace (cyclic permutation inside the trace does not 
change it) one can get 



Itt (a-Vaa^) = 

1 (cos20)A^ - (sin20)A^ 



Af + A^ 



(3.6) 



Now, we should take care about the thermodynamic limit 
(N —> 00). It is easy to see, that in this limit Eq. (|3 .6[) 
leads to (suppose that Ai is greater than A2) 



(aj) = icos2<^= -M a , 



(3.7) 



where M a is the partial magnetization corresponding to 
the subsystem with cr-spins, 

yj (e? H Z ++ - e-^Z—f + 4Z+_Z_ + 
Thus, we can express the entries of the orthogonal matrix 
diagonalising the transfer matrix in terms of the subsys- 
tem magnetization M a : 



1 f y/T+ M a -VI - M~ a 



V2 V VI - M a VI + M a 



(3.9) 



These quantities are the building blocks for calculation 
of all correlation function for the HIC. As will be clear 
from the further consideration, any correlation function 
of HIC has the form of the trace of a certain product of 
transfer matrices and special matrices representing the 
spin operators. Moreover, these matrices are constructed 
from the local expectation values of spin operators with 
the aid of the matrix A in the standard way P = A _1 PA. 
Let us write down the entries of a general matrix of that 
type 



P 



++ 
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2 

1 r 

2 

1 r 

2 
1 

2 



1 + M a ) P++ + yj\-M* (P+- + P-+) + (1 - M a ) P— 
VI - Ml (P__ - P++) + (1 + M a ) P+_ - (1 - M a ) P_+ 
VI - Ml (P__ - P++) + (1 + M a ) P_ + - (1 - M a ) P+- 
(1 - M a ) P++ - y/l - Ml (P+_ + P_+) + (1 + M CT ) 

I 



For instance, the matrix a z = A 1 cr 2 A which is the tions of aj has a particulary simple form: 
building block for constructing various correlation func- 



1 ( M a 

2 \-^Ml 



-M„ 



(3.10) 



(3.11) 
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The simple two-site correlation function is just a ther- distance r form each other: 
mal average of the product of two a spins situated at a 



(a) n=l 

I 



Thus, using the matrix form of o z and taking a thermo- 
dynamic limit one obtains: 

Wi'H-r) = J (Mi + (1 - Mi) e-i) , (3.12) 

where £ is a correlation length. So far, all results have the 
same form as those for a simple Ising chain, the difference 

I 



is only in the form of the magnetization for cr-spins, which 
in our case is given by Eqs. (|3.8[) and (|2.4[) . Now, let 
us turn to the more complicated correlation functions of 
the quantum spins with each other and mixed correlation 
functions, 



<3"X^> = £ £ St a Sl r>b J] e-^f (3.13) 

O) W n=i 

= — Tr ^T Ar - r - 1 S^T , '- 1 S^ = -Tr (A Ar - r - 1 S^A r ~ 1 Sf) , 
1 w 

(<x) (t) n=l 

= -^Tr ( T JV -'"- 1 S°T r (T 2: ) = ^Tr (a^-'-^A 1 ^) 

(a) (r) n=l 

= ^Tr (T JV - r - 1 S^T r - 1 r 6 ) = ^Tr (a^-^A^ 1 ^) , 

(<r) (t) «=1 

= Itt (t^-Vj-V,) = ^Tr (A^-^A^n) , 

(<x) (t) n=l 

= iTr (T^'-VJV) = iTr (A JV -^ 1 f a A''a) , 



where the new 2 by 2 matrices and r a are the clas- spins from the corresponding quantum clusters. Their 
sical transfer-matrix formalism representation of the a- 
component of the a-th quantum spin and classical a-th 
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entries are found according to 

fa) 

W,,, +1 =^ f, ' ,+,1+l) E T ^^ ( 3 - 14 ) 

fa) 



As the expressions under the trace in Eq. (|3.13[) have the 
same structure, let us write down a general expression 
for such a correlation function 



— Tr | 



(3.15) 



-r/C 



where the correlation length £ is the same as in Eq. (13.121) . 
Obviously, the part independent of r is the product of ex- 
pectation values of quantities corresponding to matrices 
P and R, 



(P i) = ^(A^p) = ip, 



(3.16) 



The product of two eigenvalues of the 2 by 2 transfer 
matrix is its determinant, thus, eventually one can write 
down the following general form of pair-correlation func- 
tion for HIC: 



(PjRj+r) 



—TWA 

z V 

P+-R- 



N-r-l 



PA r - 1 R (3.17) 



± e -r/t 



detT 

If one of the matrices is a the result is different, because 
the corresponding expression contains one more transfer 
matrix, 

1 



: Tr (A 

1 

2AT 



, N-r-U 



PA r ?) 

P++M a - P+_ v/l-Mje"^ 



(3.18) 



The most important conclusion one can make from this 
general consideration is the existence of only one corre- 
lation length for the any kind of HIC. Despite the inho- 
mogeneous nature of such kind of systems, the decay of 
correlations between various spins belonging to different 
blocks is characterized by the same correlation length, 
defined by the ratio of the maximal and the next to the 
maximal(in our case of 2 by 2 transfer matrix, just the 
" first" and the " second" ) eigenvalues of transfer matrix. 
This is, however, the true only for a long range correlation 
with respect to the block structure of the system; inside 
one block short range correlations are non-uniform. 
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FIG. 1: Heisenberg-Ising diamond chain. Solid (dashed) lines 
denote quantum XXZ (classical Ising) interaction bonds. All 
dashed bonds correspond to coupling constant K, the inter- 
action between quantum spins is characterized by coupling 
constant J and anisotropy A. 



several examples. As the technique, despite its clear- 
ness and simplicity, leads to rather cumbersome expres- 
sions, we will present explicit calculation for the sim- 
plest models which have been investigated earlier in a 
context of magnetic and thermodynamic properties^—. 
One can distinguish at least two classes of HIC: the mod- 
els where each block has a left-right symmetry with re- 
spect to interaction with a spins^"— i n i 14 , which implies 
that the term describing the interaction between block 
spins and their left and right cr-ncighbors has the form 
A'(<Tj +crj+i)^ a 'S| (1 , which is completely symmetric 
with respect to the permutation of the operators Sj a , and 
the models with are non symmetric with respect to left 
and right a-spins on each bloc k 10 ' 13 ' 19 . It is worth men- 
tioning, that the experimentally obtained single chain 
magnet compound with Ising and Hcisenbcrg bonds, 
which is an existing example of HIC, belongs to the sec- 
ond typ o 32 ' 34 . Magnetic properties of the asymmetric 
model are vast and complicated, they can exhibit a large 
number of ground states with a breaking of translational 
symmetry or doubling of the unit cell(block) 10 i 13 ' 19 . In 
the most symmetric case the block Hamiltonian is com- 
pletely symmetric with respect to permutation of all spins 
in it. Such models are the simplest ones. For instance, 
the various variants of the Heisenberg-Ising diamond 
chai n 4 ' 8 ' 14 and Hciscncbrg- Ising tetrahcdral^ and doublc- 
tctrahcdral chain o 9 ' 11 belong to this class of HIC. It is 
obvious, that in this case the matrices are the same 
for all spins inside one block, or, in other words, and in 
this case one does not need the additional index a. 



A. Symmetric diamond-chain 

The symmetric Heisenberg-Ising diamond chain is the 
simplest example of HIC (Sec Figure [T]). Its block Hamil- 
tonian has the form: 



IV. EXAMPLES 



H-n — </ (A (S n l S n 2 + + S n l S n 2 ) 

-(H-K (a n + a n+1 )) (S* A + S^ 2 ) (4.1) 



In this section we are going to illustrate the tech- 
nique of calculation of correlation functions for HIC on The corresponding transfer-matrix is 
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iDC 



2e 



-Pi 



e* 9 "? (ch [j3 (H — K)] + Q DC ) ch(pH) 
ch{f3H) + Q DC e - ^"? (ch [P (H + K 



DC 



Q 



DC) 



(4.2) 



where 



■)DC 



(4.3) 



As the block Hamiltonian is SO(2)-invariant with respect 
to the z-axis the corresponding E"-matrix is non-zero 
only for a = z and, as discussed above, in virtue of the 
symmetry with respect to permutation of S-operators for 
one block Xf = Sf, = X 2 . The matrix X 2 has a rather 



simple form: 



e^f sh ]fi {H 
sh (fill) 



K)} 



sh (PH) 
f sh \p (H - 



(4.4) 



K)\ 



According to Eqs. (|3.10j> - (|3.18|> . the non-zero two-site cor- 
relation functions for this model are 



2 sh (pH) M a + [e^ K ch (/3f ) - e~ f3K ch (P^f)] y/1 - M* 



16 (sh 2 (pK) + 2Q DC (ch {pK) - 1) ch [fiH)) 



-e t 



(4.5) 



i S j,a a j+r) ~ ( S j,a){ a j) = _ " 



-Pi {2sh(/3i7)A/ CT yT _ [e^ch(^f) -e-^ch(^)] (l-M 2 ) 2 } 



4Af 



DC 



e i 



r 



where 

{Sl a ) = ^[A N - 1 ^)^ l -M s (4.6) 

is the partial magnetization for the sublattice of quantum 
spins. M a is the magnetization of cr-sublattice obtained 



J 



according to Eq |3.8l from transfer-matrix (|4.2[) and Af 
is the largest eigenvalue of the transfer-matrix (|4.2|) given 
by 



Af c = e 



'* {e-^ ch (f}*f) + (2Q^ C + e/ K ) ch 



(4.7) 



4 (ch (£#) + Qocy 



r 



B. Asymmetric sawtooth chain 

The next example we are going to consider is the 
Heisenberg-Ising sawtooth chain introduced in Refs. [T3 
and [H. The lattice is depicted in Figure [5] This model 
is quite convenient for illustrating the asymmetric role of 
quantum spins from one block as well as for the demon- 
stration of the role of the additional classical spin in- 
cluded in the block. Here the block Hamiltonian contains 
two parts, one with quantum operators and another one 



- without them. Thus, one can write down the block 
Hamiltonian for the sawtooth HIC with 2 quantum spins 
in one block in the following way: 

Ul c = HI + Ka n r n - Hr n , (4.8) 

+ {H- K (<t„ + (j n+l + t„)) S* A - (i? - Ka n+ i) S^ 2 . 

Here all notations are explained in Figure [3] In con- 
trast to the case considered previously in Refs. [13 and 
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FIG. 2: Heisenberg-Ising sawtooth chain with a pair of quan- 
tum spins in one block. The unit cell consists of 4 spins, 
Sj,i, Sj,2, &j and Tj. All bonds denoted by dashed lines are 
of Ising type and the corresponding coupling constant is K, 
while the XXZ-interaction corresponding to quantum bond 
(solid line) is given by coupling and anisotropy J and A. 



Ii~3l . here, for the sake of simplicity, we consider all Ising 
bonds with the same coupling constant K. As the block 
Hamiltonian contains also the Ising spin r, in order to 
obtain the entries of transfer matrix one needs, after tak- 
ing the trace over quantum states of the pair of S-spin, 
to sum over the two possible values of the classical spin 
t. The entries of the resulting classical transfer matrix 
are: 



where 



rpSC 
1 + + 



T sc 



rpSC 

1 -+ 



2e -H jg^-S) [ch(/3 (H - K)) + Q% c ] +eP* cbY/? 



2e-Pi\e^ H -%) [ch ((3H) + Q s 2 c ] + e ~U H -f) 



ch jS [H 



K 
K 



T 



sc 



2e~i j e -£(*+*) [ch((3H) + Ql c ]+eU H +f) 
2e ~H L-f>{?+*) [ch {fi (H + K)) + Ql c ] + 



ch L3 [H - 



>f c = ch (/3 



)SC _ a pi 



ch /3 



,/A 
2 

VK 2 + J 2 A 2 



(4.10) 



The block Hamiltonian (|4.8[) is not symmetric with re- 
spect the to permutation of the quantum spin operators 



K 



r 



sc 



Qt 



} sc 



} sc 



(4.9) 



S„.i and S raj 2, thus, there are two matrices represent- 
ing the z-componcnt of the corresponding spin-operators 
using for the calculation of the correlation functions. In 
virtue of the SO(2)-invariancc of Eq. (|4.8[) the correspond- 
ing matrices for other spin components identically equal 
to zero. Let us write down the entries of the Sf and Ef 
matrices for the Heisenberg-Ising sawtooth chain. 



where 



(£?, 2 )__ 



-Pi \ eP* ch f [3 f H - |- j j + e^( H "T ) ( c h ((3 (H — K)) =p W) \ , 

~ y 

r 



4 | e #(«+f ) C h ^ ftf - ^ + e S( B +Z) (ch (/3#) ± M/) 



(4.11) 



W = K- 



i sh 



V^ 2 + J 2 A 2 



(4.12) 



and uppcr(lowcr) sings correspond to S^Sg)- The ma- 
trix t in terms of which the expectation values and cor- 
relation function of the classical spin from block Hamil- 
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tonian are expressed has the following entries: 



e -Pi \ A H ~*) [ch [fi (H - K)) + Q S 2 C ] - 



e 2 



e-H\ e U H -§) [ch((3H) + Q 
-#("+#) [ch(PH) 
?i \ -e-P( H+£ t) [ch 08 (H + K)) + Q S 2 C ] + ef* 



e 4 



ch - 
ch \B \H 4 
ch(j3[H 



K 

~2 

K 

~2 
K 
~2 



Qt 



sc 



Qt 



sc 



} sc 



(4.13) 



The all possible two-spin correlation functions (beside p,13|) ) of the system are listed below, 
the era, which is uniform for all HIC and is given by Eq. 



2 (ELi [A\ ch + Bi sh 03«)]) fe^ [At ch + B* sh 
(SUS z j+rl ) - (Sh) 2 = ± 3 L e -j 4.14 

„ _ (gLi [^cM^) +ggah W)]) (gti ^ch(^) +^sh(/3^)] / 

_ (gti [AIM^) + Blsh(f3if)]) (gtx [4ch(/3^) +^sh(/3^)]) ,. 

{bj^aj+r) - {i> jt x){crj) = ^-sc e « , 

v (?M(ffl g^VjHMgk (^ch(/3^) +^sh(/3M)) 

v , wx (ELi^ch^+^sh^jj^J^eh^+Be,^^)^ ^ 

(5? iTj+r) - (5fi)(r,-) = 5 — -e «, 

Jl1 h,m " 16£ 2 fe= oAcCh(/3fctf) 



(EL 1 [^cM^)+^^(^¥)])(Et=o[4cM/3¥)+ 5 ^M^¥)] y 

J ' J 3 ' 3 WJ2 k=0 D k ch((3kH) 

(aT v , ff w T v g^j V 71 ^? Etc Ch (gg) + g sh (gg)] 

, , N2 (Et=o[^ch(^)+55sh(^)])(Eto[4ch(^)+^ S h(^)] 

(TjT j+r ) - {TjY = - 5 — -e ( , 

3 3 3 16J2l =0 D k ch(/3kH) 



where the largest eigenvalue of the transfer-matrix (|4.9[) 
is 
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k=0 



S k ch (pkH) - 4$^ D fc ch (pkH) 



(4.15) 



\k=0 



k=0 



The coefficients in Eqs. (|4.14|) and f|4.15[) are 



At 



D 

So z 



,K 



(l + 2WM a ), A\ = 2e?- sh (p-)^l- ~M%, A3 = 0, A\ 



e-^VT^Mf,(4.16) 
P^ + e -n) Ma -2e-^W, B 1 2 =2We~^y/l^M^, B 1 3 =2e~^M <T , B\ = 0, 



:11K 



K 



B\ = (e^ +e-^)M a + 2e-^W, B\ = -2We~^ 1F%, Bj = B\, B\ = 0, 



(l-2WM a ), A\ = A\, A 2 = 0, A\ = A 



°e 3 -£ + e-^(l-2WM a ), A\=Al A 3 3 = 0, A% = A\, 



d3 r>2 r?3 nl d3 nl n3 n 



(1 + 2WM,), A| = A^, A£ = 0, A* = A 



= B-1 1 , B\ = -B\, B\ = B\, S| = 0, 



2e-^Qf c + 2e^QfC 



4° 



0, A^2e- /3 TAf (7 , 
2e ^Qf c - V / 1 3 M? , B\ = (l + 2Q S 2 C ) - 2e?$ - , 



2 (e^ ch ( (3 



,K 



A\=A\, A* = 0, A%=A\, 4 = 0, 

tp6 n5 d6 n5 p6 / >5 jp6 n d6 p5 



e""* sh 2 + (e^ - 1) Qf c - 2Qf c sh 



A" 



2e~^ (ch (£K") - 1) Ql c + ch ( p 



,3K 



ch /3 



K 



D 2 = 2Ql c sh ( P 



-PK 



i)0f c , 



+ 2e^Qf c , 5x = 2 (e"^Qf c + ch 



5 2 



r 



CONCLUSION 



In the present paper we elaborated a technique for 
calculating correlation functions in HIC. The tech- 
nique is based on the classical transfcr-matric formal- 
ism which proves itself as a most efficient and straight- 
forward method of treating the HIC, allowing one to 
obtain exact expressions for the partition function and, 
thus, for all thermodynamic functions in a quite simple 
way^ -11 ' 13-17 ' 19-22 ' 26 . However, the classical transfer- 
matrix formalism is also very convenient for the calcu- 
lation of the correlation functions of various kind, be- 
cause all statistical-mechanical expectation values are 



expressed in terms of traces of certain products of the 
transfer-matrices with matrices representing the micro- 
scopic variables in the transfer-matrix formalism^!. We 
presented a detailed scheme of calculation of the cor- 
relation function for HIC, obtained general expressions 
for the corresponding trace formulas and presented the 
general expressions for the matrix representations of the 
main microscopic variables for the typical HIC, quan- 
tum spin operators, additional classical spins from block 
Hamiltonians and intermediate spin. One of the main 
results shows that the spatial decay of any correlation 
between any pair of spin is characterized by a single cor- 
relation length obtained in a standard way as the recip- 
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rocal of the logarithm of the ratio of largest and second 
largest eigenvalues of the classical transfer matrix. In the 
simplest case considered here when intermediate classical 
cr-spins are 1/2 spins the transfer matrix has dimension 2 
by 2 and, thus, only two eigenvalues. However, the pref- 
actors are strongly affected by the microscopic variables. 
We illustrated our general formalism on two examples, 
the simplest HIC, symmetric Heiscnberg-Ising diamond- 
chain and asymmetric Heiscnberg-Ising sawtooth chain 
with a pair of quantum spins in the unit cell. We pre- 
sented the exact form of the transfer-matrix, its eigenval- 
ues, matrix representation corresponding to the classical 
transfer-matrix formalism for quantum spin operators, 
matrix representation for additional classical spin from 
block Hamiltonian and all possible pair correlation func- 
tions. A very important feature of the correlations func- 
tion of the HIC considered in the present paper is the 
absence of correlations between transverse components 
of quantum spins from different blocks. The correspond- 
ing matrix representation of the i and S% a operators 



are identically zero. This fact could be connected with 
the rotational symmetry of the block Hamiltonian with 
respect to the magnetic field axis. 

In is also worth mentioning, that the technique elab- 
orated and illustrated in the present paper, beside of 
purely theoretical interest, can be used in the practical 
calculation of the neutron diffraction spectra for the re- 
cently obtained coordination-polymer compounds which 
are the physical examples of HIG22r— 
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